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I.  i»WACT  rCanWana  an  aaatoaa  al»  H naaaaaaiy  id  Ar  Maa* 

^This  report  discusses  the  development  of  a siiqjle  numerical  model  for  the 
prediction  of  coastal  inlet  velocities,  discharge,  and  resulting  bay  level 
fluctuations.  The  model  is  a time-marching  model  that  simultaneously  solves 
the  area-averaged  momentum  equation  for  the  inlet  and  the  continuity  equation 
for  the  bay.  It  is  assumed  that  the  bay  surface  elevation  remains  horizontal 
as  it  rises  and  falls.  At  each  time  step  the  geometric  and  hydraulic  factors 
describing  the  inlet-bay  system  are  calculated  by  evaluating  flow  conditions  ^ 

(Continued) 


throughout  the  inlet  and  by  spatially  integrating  this  information  to  determine 
coefficients  of  the  first-order  differential  equations. 


This  model,  which  includes  the  important  terms  in  the  equation  of  motion, 
is  flexible,  easy  and  inexpensive  to  use,  and  gives  a good  estimate  of  the 
inlet-bay  system  hydraulics  for  various  conditions.  The  model  can  be  used  for 
single  or  multiple  inlets,  bays,  and  seas. 

This  report  includes  the  model  theory  and  derivation,  a FORTRAN  computer 
program  for  solving  the  model  equations,  and  instructions  for  use  of  the  pro- 
gram. Examples  are  given  to  illustrate  how  the  model  may  be  used  to  predict 
coastal  inlet  response  to  astronomical  tides,  seiching,  tsunamis,  and  storm 
surges.  A 


FOREWORD 


This  report  was  prepared  as  one  of  a series  of  reports  from  the 
Corps  of  Engineers'  General  Investigation  of  Tidal  Inlets  (GITI).  The 
GITI  research  program  is  under  the  technical  surveillance  of  the  Coastal 
Engineering  Research  Center  (CERC) , and  is  conducted  by  CERC,  the  U.S. 

Army  Engineer  Waterways  Experiment  Station  (WES),  other  Government  agen- 
cies, and  private  organizations.  The  model  described  in  this  report  is 
the  latest  in  a series  of  developments  of  numerical  inlet  hydraulic  models 
beginning  with  Keulegan  in  1967  and  continuing  with  recent  work  at  WES  by 
Huval  and  others. 

The  report  was  prepared  by  William  N.  Seelig,  Coastal  Structures 
Branch,  D.  Lee  Harris,  Chief,  Coastal  Oceanography  Branch,  and 
Barry  E.  Herchenroder,  Coastal  Oceanography  Branchy  CERC.  Development 
of  this  numerical  model  and  report  preparation  were  supervised  by 
R.M.  Sorensen,  Chief,  Coastal  Structures  Branch.  Technical  assistance 
was  provided  by  C.  Mason.  Civilian  members  of  the  Coastal  Engineering 
Research  Board,  Dean  M.P.  O'Brien,  Prof.  R.G.  Dean,  and  Prof.  R.L.  Wiegel, 
reviewed  this  report. 

Technical  Directors  of  CERC  and  WES  were  T.  Seville,  Jr.,  and 
F.R.  Brown,  respectively. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  166,  79th 
Congress,  approved  31  July  1945,  as  supplemented  by  Public  Law  172,  88th 
Congress,  approved  7 November  1963. 


J^  L.  CANNON 
Colonel,  Corps  of  Engineers 
Commander  and  Director 
Waterways  Experiment  Station 


^OHN  H.  COUSINS 
'Colonel,  Corps  of  Engineers 
Commander  and  Director 
Coastal  Engineering  Research  Center 
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PREFACE 


1.  TTie  Corps  of  Engiheers,  through  its  Civil  Works  program,  has  spon- 
sored, over  the  past  23  years,  research  into  the  behavior  and  characteris- 
tics of  tidal  inlets.  The  Corps'  interest  in  tidal  inlet  research  stems 
from  its  responsibilities  for  navigation,  beach  erosion  prevention  and 
control,  and  flood  control.  Tasked  with  the  creation  and  maintenance  of 
navigable  U.S.  waterways,  the  Corps  dredges  millions  of  cubic  yards  of 
material  each  year  from  tidal  inlets  that  connect  the  ocean  with  bays, 
estuaries,  and  lagoons.  Design  and  construction  of  navigation  impi'ove- 
ments  to  existing  tidal  inlets  are  an  important  part  of  the  work  of  many 
Corps'  offices.  In  some  cases,  design  and  construction  of  new  inlets  are 
required.  Development  of  information  concerning  the  hydraulic  character- 
istics of  inlets  is  important  not  only  for  navigation  and  inlet  stability, 
but  also  because  inlets,  by  allowing  for  the  ingress  of  storm  surges  and 
egress  of  flood  waters,  play  an  important  role  in  the  flushing  of  bays 
and  lagoons. 

2.  A research  program,  the  General  Investigation  of  Tidal  Inlets 
(GITI),  was  developed  to  provide  quantitative  data  for  use  in  design 

of  inlets  and  inlet  improvements.  It  is  designed  to  meet  the  following 
objectives: 

To  determine  the  effects  of  wave  action,  tidal  flow,  and  related 
phenomena  on  inlet  stability  and  on  the  hydraulic,  geometric,  and 
sedimentary  characteristics  of  tidal  inlets;  to  develop  the  knowl- 
edge necessary  to  design  effective  navigation  improvements,  new 
inlets,  and  sand  transfer  systems  at  existing  tidal  inlets;  to 
evaluate  the  water  transfer  and  flushing  capability  of  tidal 
inlets;  and  to  define  the  processes  controlling  inlet  stability. 

3.  The  GITI  is  divided  into  three  major  study  areas:  (a)  inlet 
classification,  (b).  inlet  hydraulics,  and  (c)  inlet  dynamics. 

a.  Inlet  Claaeifiaation.  The  objectives  of  the  inlet  classifi- 
cation study  are  to  classify  inlets  according  to  their  geometry,  hydrau- 
lics, and  stability,  and  to  determine  the  relationships  that  exist  among 
the  geometric  and  dynamic  characteristics  and  the  environmental  factors 
that  control  these  characteristics.  The  classification  study  keeps  the 
general  investigation  closely  related  to  real  inlets  and  produces  an 
important  inlet  data  base  useful  in  documenting  the  characteristics  of 
inlets. 


b.  Inlet  Hydvaulice.  The  objectives  of  the  inlet  hydraulics 
study  are  to  define  tide-generated  flow  regime  and  water  level  fluctua- 
tions in  the  vicinity  of  coastal  inlets  and  to  develop  techniques  for 
predicting  these  phenomena.  TIte  inlet  hydraulics  study  is  divided  into 
three  areas:  (1)  idealized  inlet  model  study,  (2)  evaluation  of  state- 
of-the-art  physical  and  numerical  models,  and  (3)  prototype  inlet 
hydraulics. 
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(1)  The  Idealized  Inlet  Model.  The  objectives  of  this 
■odel  study  are  to  determine  the  effect  of  inlet  configurations  and 
structures  on  discharge,  head  loss  and  velocity  distribution  for  a num- 
ber of  realistic  inlet  shapes  and  tide  conditions.  An  initial  set  of 
tests  in  a trapezoidal  inlet  was  conducted  between  1967  and  1970.  How- 
ever, in  order  that  subsequent  inlet  models  are  more  representative  of 
real  inlets,  a nuaber  of  "idealized"  models  representing  various  inlet 
morphological  classes  are  being  developed  and  tested.  The  effects  of 
jetties  and  wave  action  on  the  hydraulics  are  included  in  the  study. 

(2)  Evaluation  of  State-of-the-Art  Modeling  Techniques. 

The  objectives  of  this  part  of  the  inlet  hydraulics  study  are  to  deter- 
mine the  usefulness  and  reliability  of  existing  physical  and  numerical 
modeling  techniques  in  predicting  the  hydraulic  characteristics  of  inlet- 
bay  systems,  and  to  determine  whether  simple  tests,  performed  r^idly  and 
economically,  are  useful  in  the  evaluation  of  proposed  inlet  improvements. 
Masonboro  Inlet,  North  Carolina,  was  selected  as  the  prototype  inlet  whidi 
would  be  used  along  with  hydraulic  and  numerical  models  in  the  evaluation 
of  existing  techniques.  In  September  1969  a complete  set  of  hydraulic 
and  bathymetric  data  was  collected  at  Masonboro  Inlet.  Construction  of 
the  fixed-bed  physical  model  was  initiated  in  1969,  and  extensive  tests 
have  been  performed  since  then.  In  addition,  three  existing  numerical 
models  were  applied  to  predict  the  irlet's  hydraulics.  Extensive  field 
data  were  collected  at  Masonboro  Inlet  in  August  1974  for  use  in  evaluat- 
ing the  capabilities  of  the  physical  and  numerical  models. 

(3)  Prototype  Inlet  Hydraulics.  Field  studies  at  a nunber 
of  inlets  are  providing  information  on  prototype  inlet-bay  tidal  hydraulic 
relationships  and  the  effects  of  friction,  waves,  tides,  and  inlet  mor- 
phology on  these  relationships. 

a.  Inlet  Dynamioa,  Tbe  basic  objective  of  the  inlet  dynamics 
study  is  to  investigate  the  interactions  of  tidal  flow,  inlet  configura- 
tion, and  wave  action  at  tidal  inlets  as  a guide  to  improvement  of  inlet 
channels  and  nearby  shore  protection  works.  The  study  is  subdivided 
into  four  specific  areas:  (1)  model  materials  evaluation,  (2)  movable- 
bed  modeling  evaluation,  (3)  reanalysis  of  a previous  inlet  model  study, 
and  (4)  prototype  inlet  studies. 

(1)  Model  Materials  Evaluation.  This  evaluation  was  initi- 
ated in  1969  to  provide  data  on  the  response  of  movable-bed  model  mate- 
rials to  waves  and  flow  to  allow  selection  of  the  optimum  bed  materials 
for  inlet  models. 

(2)  Movable-Bed  Model  Evaluation.  The  objective  of  this 
study  is  to  evaluate  the  state-of-the-art  of  modeling  techniques,  in 
this  case  movable-bed  inlet  laodeling.  Since,  in  many  cases,  movable-bed 
modeling  is  the  only  tool  available  for  predicting  the  response  of  an 
inlet  to  improvements,  the  capabilities  and  limitations  of  these  models 
must  be  established. 


4.  This  report  describes  a numerical  model  that  can  be  used  to 
predict  inlet  channel  velocities  and  discharge  as  well  as  the  resulting 
bay  surface  level  oscillations  for  inlets  responding  to  the  tide  and 
other  long  wave  excitation.  It  has  been  developed  as  an  easy  to  use, 
inexpensive  method  for  a good  "first  look"  analysis  of  inlet  hydraulics. 
The  need  for  this  model  arose  during  preliminary  attempts  to  apply  an 
earlier  model  developed  under  the  GITI  numerical  model  evaluation 

(see  3,b,(2)  above),  and  described  in  Appendix  4 to  GITI  Report  6, 

"A  Simplified  (Lumped  Parameter)  Numerical  Simulation"  (Huval  and 
Wintergerst,  1977).  Although  both  models  expand  upon  the  Keulegan  (1967) 
concept  of  a simple  one-dimensional  analysis,  this  model  supersedes  that 
of  Appendix  4 since  it  provides  more  accurate  results  and  is  applicable 
to  a wider  variety  of  situations  and  conditions. 

5.  Included  in  this  report  are  a derivation  of  the  numerical  model, 
dociimentation  of  the  FORTRAN  computer  program  used  to  apply  the  nx>del, 
and  example  applications  of  the  model  to  evaluate  the  hydraulic  conditions 
at  selected  tidal  and  nontidal  inlets. 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUIUiMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.  39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

mi  les 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  X 10" 

^ kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins^ 

^To  obtain  Celsius 

(C)  tciiipcraturc 

readings  from  Fahrenheit  (F)  readings. 

use  formula:  C = 

(5/9)  (F  -32). 

To  obtain  Kelvin 

(K)  readings,  use 

formula:  K = (5/9)  (F  -32)  + 273.15. 

i 
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SYMBOLS  AND  DEFINITIONS 


Symbol 
in  the 
report 

Symbol  in  the 
computer  program 
INLET 

Definition 

A(I.J) 

cross-sectional  area  (square  feet) 

Ai 

AB 

cross-sectional  area  at  the  bay  end  of 
the  inlet  (square  feet) 

Ag 

AS 

cross-sectional  area  at  the  sea  end  of 
the  inlet  (square  feet) 

^bay 

ABAY 

bay  surface  area  (square  feet) 

Ac 

ABY 

bay  surface  area  at  datum  (square  feet) 

- 

AO 

ocean  forcing  amplitude  (feet) 

B 

B(I,J) 

width  (feet) 

Ci 

C(I) 

flow  resistance  parameter 

c,.c, 

C1,C2 

coefficients  to  evaluate  Manning's  bottom- 
friction  factor,  n,  where  n = Cj  - C2  * 

D 

D(I.J) 

total  water  depth  (feet) 

^bay 

— 

depth  of  the  bay  (feet) 

‘Wix 

— 

maximum  water  depth  in  the  inlet  (feet) 

F 

F 

total  inlet  friction 

— 

minimum  friction  function 

Fx 

— 

friction  in  the  grid  cells  of  a cross 
section 

g 

G 

acceleration  of  gravity  (feet  per  second 
squared) 

h 

H(I.J) 

water  level  above  datum  (feet) 

HB 

water  level  at  the  bay  end  (feet) 

HS 

water  level  at  the  sea  end  (feet) 

HOURS 

model  time  (hours) 

SYMBOLS  AND  DEFINITIONS- -Continued 


IC 

IC 

number  of  channels  in  a grid 

— 

geometry  integral  (see  text) 

IS 

IS 

number  of  sections  in  a grid 

I.J 

subscripts  indicating  grid  location 

k 

— 

a constant  in  Manning's  equation  to 
maintain  consistent  units 

^bay 

— 

length  of  a bay  (feet) 

L(I.J) 

length  of  a grid  (feet) 

^in 

LENGTH 

inlet  length  (feet) 

M 

N INLETS 

number  of  inlets  connecting  the  bay  to 
the  sea 

m 

NI 

subscript  identifying  the  inlet  number 

n 

N(I.J) 

Manning's  bottom  friction  coefficient 

- 

NCYCLES 

number  of  forcing  water  level  cycles  used 
in  computation 

— 

NT 

number  of  time  steps  used  in  computations 

Q 

Y(NI)  p 

discharge  of  an  inlet  (cubic  feet  per 
second) 

Qij 

Q(I.J) 

discharge  of  a grid  (cubic  feet  per 
second] 

^nflou 

QINFLO 

net  discharge  of  water  into  the  bay  from 
sources  other  than  the  inlet  (cubic  feet 
per  second) 

Qj< 

QT 

total  discharge  of  all  inlets  connecting 
a sea  to  a bay  (cubic  feet  per  second) 

— 

QINT(NI) 

estimated  inlet  discharge  at  time  the 
model  starts  (cubic  feet  per  second) 

R 

— 

remainder  terms  (neglected) 

Tf 

T 

forcing  wave  period  in  the  sea  (hours 
or  seconds) 

II 


SYMBOLS  AND  DEFINITIONS— Continued 


THEIX 

X 

THTF 

TIME 

V(I,J) 

VBAR 

W(IJ,) 

DELT 

BETA 

ZETA 

DERBY 


inlet-hay  Helmholtz  period  (hours  or 
seconds) 

inlet-bay  Helmholtz  period  estimated  by 
neglecting  inlet  friction  (hours  or 
seconds) 

time  of  model  operation  (seconds) 

ratio  of  the  Helmholtz  to  forcing  wave 
period 

tine  (seconds) 

water  velocity  (feet  per  second) 

cross-sectional  mean  water  velocity 
(feet  per  second) 

mean  water  velocity  across  the  minimum 
area  section  (feet  per  second) 

grid  weighting  function  for  distributing 
flow  throughout  an  inlet  grid  flow  net 

distance  along  a channel  (feet) 

bay  limit 

sea  limit 

distances  perpendicular  to  the  main  axis 
of  the  inlet  channel  (feet) 

variable  finite-difference  time  step 
(seconds) 

bay  surface  area  variation  parameter 
relating  bay  area  to  bay  water  level 
where:  ABAY  - ABY  (1.  ♦ BETA  * HB) 

inlet  side  slope,  dZ/dy,  where  Z is 
elevation 

Lagrangian  multiplier 

derivative  of  inlet  discharge  with  respect 
to  time  (cubic  feet  per  second  squared) 

component  of  the  stress  tensor 
representing  the  bottom  stress 


A SPATIALLY  INTEGRATED  NUMERICAL  MODEL  OF  INLET  HYDRAULICS 


by 

william  N.  Seelig,  D.  Lee  Harriet 
and  Barry  E.  Herohevroder 

I.  INTRODUCTION 

Quick,  inexpensive  estimates  of  inlet  velocities  and  bay  water 
surface  levels  for  tidal  or  nontidal  sea  level  fluctuations  are  needed 
in  planning  the  design,  construction,  and  maintenance  of  coastal  inlets. 
Field  data  are  often  unavailable,  and  available  data  are  often  incomplete. 
In  addition,  hydraulic  characteristics  for  proposed  inlets  are  unavailable 
and  must  be  predicted. 

This  study  discusses  these  needs  by  developing  a numerical  model 
that  can  be  used  to  estimate  inlet  velocities,  discharge,  and  bay  water 
levels  as  functions  of  time  for  a given  time -dependent  sea  level  fluctu- 
ation. The  objective  of  this  research  was  to  develop  a model  that  could 
he  used  to  quantitatively  predict  hydraulics  for  as  wide  a range  of 
conditions  as  possible,  while  being  easy  and  inexpensive  to  use.  Inlet 
hydraulics  are  predicted  in  this  model  by  marching  through  time,  simul- 
taneously solving  the  momentum  equation  for  flow  in  the  inlet  and  the 
continuity  equation  relating  the  bay  level  to  inlet  discharge.  The  momen- 
tum equation  is  evaluated  at  each  time  step  by  integrating  two-dimensional 
information  into  coefficients  of  the  equation,  using  a weighting  function 
(see  App.  A)  and  a flow  net  (see  App.  B)  which  systematically  distribute 
flow  throughout  the  inlet. 

The  advantages  of  this  model  are  that  it  requires  a minimum  amount 
of  input  data,  and  it  is  easy  and  inexpensive  to  use  as  compared  to  a 
full  two-dimensional  model.  Other  features  of  the  model  are  that  it  in- 
cludes all  potentially  important  terms  developed  from  the  three-dimensional 
momentum  equations,  and  it  allows  for  special  situations.  For  example, 
water  level  fluctuations  in  the  sea  can  be  any  function  of  time,  the  area 
of  the  bay  can  be  a function  of  water  level,  and  inflow  into  the  bay 
from  sources  other  than  the  inlet  can  be  a function  of  time  or  water 
level  of  the  bay.  Also,  the  cross-sectional  area  of  the  inlet  can  be 
taken  as  a function  of  the  local  water  depth.  This  model  can  be  used  to 
predict  hydraulics  for  tidal  or  nontidal,  single  and  multiple  inlet 
systems.  It  assumes  that  the  bay  level  "puiq>s''  (i.e.,  rises  and  falls) 
at  the  sane  rate  and  phase  throughout  the  bay. 

II.  INLET  HYDRAULICS 

An  inlet-bay  system  typically  consists  of  a "sma"  (e.g.,  ocean  or 
lake)  connected  to  a "bay"  by  one  or  more  inlets  (Fig.  1).  Long  waves 
in  the  sea  (i.e.,  astronomical  tides,  seiches,  storm  surges,  tsunamis, 
or  other  water  level  fluctuations)  generate  the  primary  hydraulic  response 
in  the  inlet-bay  system.  The  difference  in  water  level  between  the  bay 
and  sea,  caused  by  the  sea  forcing  fluctuations,  results  in  reversing 
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currents  in  the  inlet  which  in  turn  cause  the  bay  level  to  change.  The 
bay  level  may  also  be  influenced  by  inflow  from  and  outflow  to  other 
sources  (i.e.,  river  discharge  or  evaporation]  which  may  be  a function 
of  time  or  of  the  bay  level.  Additional  factors  such  as  wind  stress, 
water  density  variations,  and  the  earth's  rotation  may  also  affect  the 
system. 

The  prediction  of  reversing  currents  in  the  inlet  and  bay  level 
fluctuations  requires  a knowledge  of  the  geometries  of  the  inlets  and 
bay,  the  water  level  fluctuations  in  the  sea  which  force  the  inlet-bay 
system,  and  a model  (analytical,  numerical,  or  physical)  to  predict  the 
system  response.  For  most  inlets  and  bays,  the  geometries  can  be  measured 
in  the  field,  obtained  from  dredging  records,  or  evaluated  from  hydro- 
graphic  charts.  Water  level  fluctuations  in  the  sea  must  be  measured 
in  the  field,  predicted  from  tide  tables  or  physical  or  numerical  hydro- 
dynamic  models. 

The  complete  three-dimensional  equations  of  motion  that  describe  the 
hydrodynamics  of  the  inlet-bay  system  (Fig.  1)  are  complicated.  Rather 
than  try  to  solve  the  complete  equations,  various  investigators  have 
simplified  them  and  attempted  to  model  only  the  factors  that  the  investi- 
gators considered  important  in  inlet  hydraulics.  Some  models,  such  as 
physical  distorted  scale  models  and  some  two-dimensional  vertically  inte- 
grated numerical  models,  have  been  used  to  adequately  predict  inlet 
hydraulics  for  most  engineering  purposes  (Harris  and  Bodine,  1977). 
However,  operation  of  these  models  is  relatively  expensive,  requires 
expert  operators,  and  usually  takes  several  months  of  effort  to  obtain 
the  required  results. 


Simple  analytical  models  such  as  those  of  Brown  (1928)  and  Keulegan 
(1967),  which  lump  important  parameters  into  a single  quantity,  require 
a number  of  restrictive  assumptions.  O'Brien  and  Clark  (1974)  concluded 
"The  Keulegan  approach  and  similar  analyses  of  inlet  hydraulics  provide 
a useful  qualitative  framework  for  ordering  data  but  they  apply  quanti- 
tatively only  to  small  inlets  and  lagoons  with  simple  inlet  channel 
geometry. " 

More  advanced  lumped  parameter  models,  such  as  those  of  King  (1974) 
and  Huval  and  Wintergerst  (1977),  have  included  additional  terms  in  the 
equation  of  motion  and  have  considered  other  effects  to  generalize  the 
models.  However,  many  of  these  more  advanced  lumped  parameter  models 
were  based  on  earlier  models  and  in  the  processes  of  expansion  some 
generality  was  lost.  Huval  and  Wintergerst  (1977)  and  Harris  and  Bodine 
(1977)  discuss  previous  models. 


III.  THE  NUMERICAL  MODEL 


1.  Assumptions 


This  spatially  integrated  numerical  model  is  based  on  the  equations 
derived  from  the  complete  equations  of  motion  by  Harris  and  Bodine. (1977) 
Assumptions  made  in  the  process  of  model  derivation  were: 


(a)  Sea  level  is  a specified  function  of  time. 
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(b)  The  bay  water  level  remains  horizontal.  This  means 
that  the  bay  water  level  rises  and  falls  at  the  same  rate  through- 
out the  entire  bay  at  each  point  in  time.  This  occurs  when  the 
length  of  the  long  wave  forcing  the  system  is  much  longer  in  the  bay 
than  the  longest  axis  of  the  bay. 

(c)  The  bay  is  connected  to  the  sea  by  one  or  more  inlets. 

(d)  At  least  one  inlet  must  continuously  connect  the  bay  to 
the  sea.  Some  areas  of  inlets  may  go  dry  during  the  water  level 
cycle,  and  one  or  more  inlets  may  go  dry  as  long  as  one  inlet 
contains  water. 

(e)  Bay  water  surface  area  is  a function  of  bay  water  level 
(or  a function  of  time). 

(f)  Inlet  cross-sectional  area  is  a function  of  local  depth 
(or  a function  of  time). 

(g)  The  local  water  level  slope  in  the  inlet  is  assumed  to  be 
linearly  related  to  the  local  friction  loss  along  the  inlet  between 
the  sea  and  bay  levels. 

(h)  There  is  a water  level  drop  along  the  inlet  that  is  pro- 
portional to  the  unrecovered  velocity  head  lost  through  turbulent 
eddy  diffusion  in  the  bay  (floodflow)  or  sea  (ebb  flow).  (Alter- 
nate schemes  could  be  used.) 

(i)  Storage  of  water  ih  the  inlet  is  negligible.  This  means 
that  the  flow  into  the  inlet  is  equal  to  the  flow  exiting  the  inlet 
at  any  time.  In  addition,  the  volume  of  water  stored  in  the  inlet 
between  high  and  low  water  should  be  small  compared  to  the  tidal 
prism.  This  is  generally  the  case  if  the  surface  area  of  the  bay 
is  much  larger  than  the  surface  area  of  the  inlet. 

(j)  Wind  stress  on  the  inlet  and  bay  surfaces  is  negligible. 
This  means  that  the  model  is  most  useful  for  cases  when  the  wind  is 
light,  has  a short  duration,  or  has  a short  fetch  over  the  bay. 

(k)  Water  has  constant  properties  throughout  the  inlet  and 
bay.  No  attempt  has  been  BUde  to  model  smltwater  intrusion  or 
other  density  gradient  effects. 

(l)  Radiation  stress  (the  interaction  with  wind  waves)  is 
neglected. 

(m)  Coriolis  effects  are  neglected. 

If  the  characteristics  of  the  inlet-bay  systems  modeled  are  lot 
consistent  with  these  assu^>tions,  model  results  may  be  in  error. 
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2.  Derivation. 


The  derivation  of  this  model  begins  with  the  one-dimensional  equation 
of  motion  as  derived  by  Harris  and  Bodine  (1977): 


where 


cross-sectional  mean  water  velocity  in  the  inlet 
(positive  on  floodflow) 


distance  along  the  main  axis  of  the  inlet 


* water  level  above  some  datum 


acceleration  due  to  gravity 


inlet  cross-sectional  flow  area  at  x 


” component  of  the  stress  tensor  at  the  bottom  of  the 
^ inlet  in  the  direction  of  the  main  axis  of  the  inlet 

R ■ remainder  terms,  which  are  neglected  in  this  model  (R“0). 

This  means  that  the  water  level  is  taken  as  constant  at 
each  inlet  cross  section.  See  Harris  and  Bodine  (1977) 
for  a discussion  of  these  neglected  terms. 


The  first  term  on  the  left  of  equation  (1)  is  the  temporal  acceler- 
ation, the  second  term  is  the  convective  or  advective  acceleration,  the 
third  term  is  the  slope  of  water  surface  along  the  inlet,  and  the  fourth 
term  is  the  bottom  stress. 

To  obtain  a simplified  equation,  the  expressions  in  equation  (1) 
are  integrated  over  the  length  of  the  inlet  between  the  sea  and  the  bay, 
where  Xm  and  x.  are  the  respective  limits: 


Carrying  out  some  of  the  integrations  and  rearranging,  equation 
(2)  becomes 


In  equation  (3),  terms  involving  3x,/3t  and  3Xg/3t  have  been 
set  to  zero  since  and  Xg  are  taken  to  be  independent  of  time. 

From  continuity  the  cross-sectional  mean  inlet  water  velocity  is 
equal  to  the  inlet  discharge,  Q,  divided  by  the  inlet  cross-sectional 
area , : 


Substituting  equation  (4)  and  using  the  product  rule  for  integration 
the  first  term  on  the  left  of  equation  (3)  can  be  integrated  to  yield: 


where  the  second  part  of  the  equation. 


3t 


is  taken  as  zero  because  channel  storage  terms  are  neglected. 

After  substitution  of  equation  (4) , the  second  and  third  terms  on 
the  left  of  equation  (3)  are 


1^ 
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(6) 


where  A,  and  Ag  are  the  cross-sectional  areas  of  the  inlet  at  the 
bay  and  sea  ends  of  the  inlet,  and  h^  and  hg  are  the  water  levels 
in  the  bay  and  sea  referenced  to  a common  datum.  The  convective  acceler- 
ation, the  left  term  of  equation  (6),  may  also  be  expressed  in  terms  of 
empirical  loss  coefficients  (see  App.  C). 

The  bottom  stress  is  evaluated  by  using  Manning's  equation; 


(7) 


where  u is  the  water  velocity  in  the  inlet,  D is  the  water  depth  at 
that  point,  and  k is  a conversion  factor  to  adapt  Manning's  equation  to 
the  system  of  units  used.  The  value  of  k is  [1.0  meter  (3.2808  feet) 
per  second  cubed] The  absolute  value  function  of  u accounts  for 
the  alternating  direction  of  bottom  stress.  The  bottom  stress  is  approxi- 
mated by  determining  the  water  velocity,  u,  at  a number  of  grid  points 
throughout  the  inlet. 

The  value  of  u is  obtained  from  a plan  view  flow  net  system  of 
channels  and  cross  sections  which  are  representative  of  flow  throughout 
the  tidal  cycle  (see  App.  B).  The  main  axis  of  each  channel  is  drawn 
approximately  parallel  to  the  direction  of  flow  and  the  cross  sections 
are  drawn  approximately  perpendicular  to  the  flow.  The  subscript,  i, 
then  denotes  the  cross-section  number  and  the  subscript,  j,  indicates 
the  channel  number  (Fig.  2).  Section  i»l  is  the  seaward  cross  section 
and  section  i*IS  is  the  inlet  cross  section  at  the  entrance  to  the  bay. 
Channel  j»l  and  channel  j»IC  are  the  two  channels  adjacent  to  land. 

A typical  cell,  denoted  as  cell  (i,j),  consists  of  that  part  of  channel 
j situated  between  cross-sections  i and  i+1.  The  water  velocity  in 
cell  (i,j),  u£j,  is  assumed  to  be  at  the  centroid  of  cell  (i,j)  and 
to  act  parallel  to  the  axis  of  channel  j (Fig.  2). 
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Figure  2.  Inlet  grid  system. 


A weighting  function,  then  used  to  determine  what  fraction 

of  the  total  inlet  flow,  Q,  passes  through  a grid  cell  (i,j)  at  an 
instant  in  time: 


where  is  the  discharge  in  cell  (i,j).  Various  weighting  functions 

are  discussed  later  in  this  section. 

The  mean  water  velocity  in  cell  (i,j),  taken  equal  to  the 

discharge  in  the  cell  divided  by  the  mean  cross-sectional  area  of  the 
cell  perpendicular  to  flow,  i.e*. 


Substituting  equations  (7)  and  (9)  into  the  fourth  term  of  equation  (3) 
and  integrating  over  x and  y give  the  total  bottom  friction,  F: 


1 


I Wjj  QlWy  Q Hj 


(10) 


where  n-in  is  the  Manning's  coefficient  of  friction  for  each  grid  cell, 
is  the  mean  instantaneous  water  depth  in  a cell,  and  and 

are  the  mean  cell  width  and  length,  respectively.  All  of  the  cell  param- 
eters are  taken  to  apply  to  the  centroid  of  the  cell  (Fig.  2). 

Let  Ig  be  a geometry  integral  defined  for  convenience  as 


which  has  units  of  length. 

Substituting  equations  (5),  (6),  and  (10)  into  equation  (3),  multi- 
plying by  Ig,  and  setting  dQ/dt  equal  to  the  other  terms  will  obtain: 
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where  is  the  instantaneous  surface  area  of  the  hay. 

There  are  several  methods  available  for  solving  the  simultaneous 
differential  equations  (eq.  12  for  M inlets,  and  eq.  14  for  a total 
of  M+1  equations).  The  method  selected  for  this  model  is  a fourth- 
order  Runge-Kutta-Gill  technique.  Advantages  of  this  method  are  that  it 
is  self-starting,  extremely  stable,  may  use  a longtime  step,  has  wide 
application,  and  converges  quickly.  The  main  disadvantage  of  this  tech- 
nique is  that  it  may  cost  approximately  twice  as  much  as  some  methods 
because  the  Runge-Kutta-Gill  technique  uses  two  calculations  for  each 
time  step  to  check  error  bounds  and  establish  the  time  step  (International 
Business  Machine,  1970).  On  a CDC  6600  computer  the  total  computer  cost 
for  computations  per  inlet  for  a tidal  cycle  using  the  Runge-Kutta-Gill 
method  has  been,  at  most,  several  dollars  (less  than  $2  per  tide  cycle) 
for  most  test  inlets.  However,  the  ease  in  use  of  this  method  justifies 
its  cost. 

The  computer  program  based  on  this  model  (INLET)  is  presented  in 
Appendix  D. 


3.  Weighting  Functions. 


The  weighting  function,  (the  fraction  of  the  total  flow  that 

passes  through  a grid  cell  at  a time  step),  provides  a systematic  method 
of  distributing  flow  throughout  an  inlet  for  use  in  evaluating  the  bottom 
stress  (see  App.  A).  There  are  three  weighting  functions  developed  in 
this  report. 

First,  a weighting  function  to  distribute  flow  between  channels  at 
each  cross  section  so  that  total  friction  in  the  section  is  minimized 
(option  IWT=1  in  the  computer  program  INLET).  This  function  may  allow 
some  water  to  move  perpendicular  to  the  main  axis  of  the  flow  net  at 
each  cross  section,  but  the  flow  should  be  small  for  a well-drawn  net. 

Second,  a weighting  function  is  developed  by  assuming  that  all  flow 
is  parallel  to  the  streamlines  of  the  flow  net  and  distributes  the  dis- 
charge in  each  channel  of  the  inlet  to  minimize  overall  friction 
(option  IWT=2  in  the  computer  program  INLET).  This  method  is  consistent 
with  the  equations  of  motion  used  in  the  model  derivation  and  assumes 
that  the  grid  system  permits  an  accurate  representation  of  the  inlet 
streamline  patterns. 

In  practice  both  of  the  minimum  friction  weighting  functions  (IWT=1 
and  IWT=2) produce  similar  results  for  many  inlets.  Comparisons  with 
prototype  measurements  show  that  either  minimum  friction  weighting  func- 
tion adequately  predicts  the  flow  distribution  across  inlet  cross  sections 
tested.  For  example,  at  Brown  Cedar  Cut,  Texas,  minimum  friction  weight- 
ing predicted  the  fraction  of  flow  at  any  point  at  a cross  section  to 
within  several  percent  of  the  measured  total  flow  (Fig.  3).  At  times  in 
the  tidal  cycle,  minimum  friction  weighting  slightly  over  or  under  pre- 
dicts flow  in  parts  of  the  section,  but  on  the  average  this  difference 
will  have  little  effect  on  the  final  result. 

Differences  between  prototype  and  minimum  friction  weighting  distri- 
butions may  occur  for  several  reasons.  First,  Manning's  uniform  flow 
friction  relation,  as  used  in  this  model,  may  not  completely  describe 
friction  losses  for  unsteady  inlet  flow.  For  example.  Manning's  equation 
may  overestimate  friction  when  discharge  is  increasing  because  turbulence 
has  not  developed  to  the  point  that  would  be  reached  for  steady  flow  of 
the  same  discharge.  As  the  magnitude  of  discharge  decreases,  the  opposite 
effect  occurs.  Turbulence  is  higher  than  for  the  same  steady-state 
discharge,  so  Manning's  equation  underestimates  actual  friction.  Changes 
in  bed  forms  throughout  the  tide  cycle  may  also  change  frictional  resist- 
ance. 

Jet  formation  in  flow  exiting  the  inlet  throat,  not  accounted  for  in 
these  weighting  functions,  may  also  cause  differences  between  observed 
and  predicted  flow  patterns.  However,  an  empirical  coefficient  may  be 
used  to  account  for  losses  caused  by  jets  (discussed  in  App.  C) . 

See  Appendix  B for  a detailed  discussion  of  drawing  flow  nets  based 
on  the  minimum  friction  weighting  function  (IWT^l). 
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The  third  weighting  function  (1WT»3)  is  the  same  as  that  used  in 
lumped  parameter  models  (Huval  and  Kintergerst,  1977).  It  assumes  that 
discharge  is  the  same  in  all  grid  cells  throughout  the  inlet  at  any  in- 
stant in  time.  This  means  that  the  flow  in  each  grid  is  equal  to  the 
total  flow  for  the  inlet  divided  by  the  number  of  channels,  IC.  The 
weighting  function,  is  evaluated  by: 

yiij  = 1.0/ IC  (15) 

Generally,  this  function  will  not  accurately  reproduce  prototype  condi- 
tions because  it  is  almost  impossible  to  visually  draw  the  grid  that 
corresponds  to  this  condition.  In  addition,  the  function  does  not 
consider  changes  in  flow  distribution  that  occur  during  the  tidal  cycle. 
However,  this  weighting  function  is  useful  in  obtaining  an  upper  limit 
for  frictional  effects  because  it  overemphasizes  friction  in  most  cases. 
By  producing  high  friction  the  function  tends  to  give  a lower  limit  for 
bay  level  fluctuations,  mean  inlet  velocity,  and  inlet  discharge. 

Weighting  functions  are  organized  as  subroutines  in  the  computer 
program  INLET,  so  that  other  weighting  functions  may  easily  be  added  to 
this  numerical  model  as  they  are  developed. 

IV.  USE  OF  THE  COMPUTER  PROGRAM 
1.  General  Conditions  for  Model  Use. 


The  computer  program  INLET  (App.  D)  estimates  inlet  hydraulics  by 
solving  equations  (12),  (13),  and  (14)  with  the  use  of  appropriate  geo- 
metric and  hydraulic  inputs.  This  section  summarizes  the  conditions 
that  can  be  modeled  consistent  with  assumptions  discussed  in  Section  III. 
Details  of  the  application  of  the  computer  program  are  discussed  in 
Appendix  D. 

The  conditions  required  when  using  the  model  are: 

(a)  The  sea  is  much  larger  than  the  inlets  and  bay. 

(b)  The  bay  is  much  larger  than  the  inlets. 

(c)  The  bay  and  at  least  one  of  the  inlets  must  contain 

water  throughout  the  water  level  cycle. 

(d)  The  forcing  seawater  level  fluctuations  must  be  speci- 
fied as  a function  of  time.  Start  the  model  when  the  sea  level 
is  zero. 

(e)  The  inflow  or  outflow  from  sources  other  than  the  inlet 
must  be  a specified  function  of  time,  a function  of  bay  water 
level,  or  constant  (the  coBq;>uter  model  is  programed  for  a 
constant  inflow) . 
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with  high  temporal  acceleration,  such  as  inlets  on  the 
Great  Lakes. 

(m)  The  program  INLET  is  designed  so  that  the  sea 
forcing  is  sinusoidal  or  is  described  by  water  levels  sampled 

uniform  rate,  ^inflow  constant,  ^bay  ® linear 
function  of  the  bay  water  level,  and  the  inlev  side  slope  is 
linearly  related  to  the  local  water  depth.  For  more  general 
conditions,  appropriate  programing  changes  should  be  made  to 
the  computer  program.  For  example,  if  the  surface  area  of  the 
bay  is  a nonlinear  function  of  bay  water  level,  the  program 
statement  that  evaluates  bay  surface  area  should  be  changed 
to  incorporate  the  given  function. 

(n)  The  recommended  method  of  using  this  model  is 
similar  to  that  for  other  inlet  models.  Obtain  prototype 
geometries,  sea  and  bay  water  levels,  and  inlet  velocities 
for  the  system.  Using  this  prototype  information,  calibrate 
the  model  so  that  the  predicted  bay  levels  and  inlet  velocities 
accurately  reproduce  the  prototype  data.  The  suggested  method 
of  calibration  is  to  vary  the  Manning's  bottom  friction  factor, 
n,  or  ebb  and  flood  entrance  and  exit-loss  coefficients 

(App.  C;)  until  good  agreement  is  obtained  between  observed  and 
predicted  hydraulics.  For  short  records  of  field  data  it  is 
recommended  that  the  model  first  be  calibrated  so  that  predicted 
inlet  velocities  or  discharges  reproduce  prototype  conditions, 
because  velocities  are  most  sensitive  to  changes  in  model  param- 
eters. Bay  levels  are  a form  of  integrated  inlet  discharge, 
so  levels  are  less  sensitive  to  variation  in  model  parameters. 

After  the  model  is  calibrated  for  velocities,  check  bay  level 
predictions. 

If  additional  prototype  data  are  available,  the  calibrated 
model  should  be  run  with  the  additional  data  to  verify  that 
the  model  produces  adequate  results. 

If  no  prototype  data  are  available  for  calibration,  use 
values  of  n and  loss  coefficients  previously  calibrated  on 
the  model  for  similar  inlets.  The  examples  presented  later  in 
this  section  provide  estimates  of  values  that  apply  to  various 
types  of  inlets.  Specific  applications  of  the  model  are  also 
presented, 

A preliminary  estimate  of  n can  be  obtained  from  the  linear 
relation  recommended  by  Masch,  Brandes,  and  Reagan  (1977)  for 
grid  cells  of  tidal  inlet  models.  They  reported  that  n.  is  weakly 
dependent  on  water  depth  by  the  relation: 

n-Cj-C^D  (18) 
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where  D is  the  Stillwater  depth,  because  the  probability  of 
vegetation  decreases  with  depth.  For  water  depths  less  than  9.1 
meters  (30  feet)  and  greater  than  1.2  meters  (4  feet)  they  recommend 
Cj  > 0.0377  and  C,  « 0.000667.  This  condition  is  assumed  if  a 
relation  for  n is  not  specified  in  input  data.  For  depths  less 
than  1.2  meters  they  recommend  “ 0.055  and  C2  ® 0.005. 

(0)  Note  that  asymmetry  of  the  inlet  and  bay  geometries 
throughout  the  forcing  cycle  will  result  in  asymmetrical  inlet 
hydraulics.  For  example,  a sinusoidal  ocean  tide  may  produce 
a nonsinusoidal  bay  tide  due  to  effects  such  as  the  change  in 
inlet  cross-sectional  area  with  water  depth  (Keulegan,  1967). 

2.  Examples. 

The  model  was  applied  to  a variety  of  tidal  and  nontidal  inlets. 

Five  examples  were  selected  to  illustrate  the  range  of  conditions  that 
can  be  modeled:  Pentwater  Inlet,  Michigan,  to  show  the  response  of  a 
simple  geometry,  nontidal  inlet  to  forcing  at  different  wave  periods 
due  to  seiching  of  Lake  Michigan;  a hypothetical  harbor,  which  illustrates 
the  application  of  the  model  to  predict  tsunami-induced  hydraulics  at  a 
harbor;  Masonboro  Inlet,  North  Carolina,  to  illustrate  the  evaluation  of 
tidal  inlet  hydraulics;  Indian  River,  Delaware,  to  show  how  the  model  can 
predict  the  effect  of  storm  surge  at  a tidal  inlet;  and  Cabin  Point  Creek, 
Virginia,  to  show  the  effect  of  adding  a second  inlet  to  a one-inlet 
tidal  system. 

a.  Pentwater  Michigan.  Pentwater  Inlet,  located  on  the 

east  coast  of  Lake  Michigan,  is  an  example  of  a nontidal  Great  Lakes 
inlet  (Fig.  4).  Like  many  Great  Lakes  inlets,  Pentwater  is  controlled 
by  parallel  jetties  and  retaining  walls,  and  is  periodically  dredged  to 
maintain  the  channel. 

(1)  Geometry.  The  channel  is  610  meters  (2,000  feet)  long  with 

a width  of  44  meters  (145  feet)  and  a minimum  depth  of  3.6  meters  (12  feet). 
A survey  of  the  channel  at  six  cross  sections  with  a spacing  of  122  meters 
(400  feet)  was  used  to  calculate  the  cross-sectional  flow  area  along  the 
inlet  (Fig.  5).  Since  the  inlet  cross-sectional  area  and  depth  are  approx- 
imately uniform,  the  inlet  was  modeled  using  a one-channel  flow  net.  The 
area  of  the  bay,  Pentwater  Lake  (1.68  w 10**  square  meters,  1.812  x 10^ 
square  feet),  was  measured  from  Lake  Survey  Chart  77. 

Figure  6 summarizes  geometric  measurements  used  in  the  Pentwater 
Inlet  model. 

(2)  Forcing  Hydraulics.  Spectral  analysis  of  water  levels  in 
Pentwater  Lake  Indicates  that  astronomical  tides  are  not  important  and 
that  Pentwater  is  being  forced  by  some  higher  seiching  modes  of  Lake 
Michigan  (Fig.  7).  For  example,  the  5. 3-hour  peak  in  the  spectra  corre- 
sponds to  the  second  longitudinal  mode  of  oscillation  of  Lake  Michigan; 
the  3.5-hour  peak  is  the  third  mode,  etc.  (Seelig  and  Sorensen,  1977). 


Spectral  Density 


Figure  7.  Spectral  analysis  of  Pentwater  Lake  water  levels 
(3,  4,  and  5 November  1974). 
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(3)  Calibration.  The  Pentwater  Inlet  model  was  calibrated  by 
using  the  geometries  in  Figure  6,  observed  water  levels  in  Lake  Michigan 
(Duane  and  Saylor,  1967)  to  force  the  model,  and  the  Masch,  Brandes,  and 
Reagan  (1977)  method  of  determining  Manning's  n.  Various  values  of  the 
flood-  and  ebb-loss  coefficients  were  tested  (see  App.  C) ; a value  of  1.0 
was  found  to  give  the  best  prediction  of  measured  bay  water  levels  (Duane 
and  Saylor,  1967).  A comparison  of  the  calibrated  results  for  a selected 
period  of  record  (Fig.  8)  shows  that  the  forcing  water  levels  contain  a 
number  of  wave  periods,  but  waves  shorter  than  the  Helmholtz  period  (the 
natural  pumping  period)  of  Pentwater,  1.8  hours,  are  filtered  out  inside 
the  harbor. 


V This  calibration  is  checked  by  comparing  results  from  sinusoidal  long 
wave  forcing  of  the  model  to  spectral  characteristics  of  level  fluctua- 
tions measured  simultaneously  in  the  bay  and  sea.  The  numerical  model 
was  run  for  four  cycles  in  these  tests  so  that  transient  terms  due  to 
startup  conditions  had  died  out.  Separate  runs  were  made  at  typical 
seiche  periods  for  Lake  Michigan,  using  a forcing  wave  amplitude  of  3 
centimeters  (0.1  foot).  The  results  for  Pentwater,  forced  at  a period  of 
2 hours,  are  shown  in  Figure  9,  which  shows  the  given  sea  forcing  levels, 
predicted  bay  levels,  and  mean  inlet  velocity  at  the  minimum  area  section. 
The  relative  importance  of  each  term  in  the  equation  of  motion,  as  a 
function  of  phase  of  the  oscillation  is  shown  in  Figure  9.  In  order  to 
emphasize  the  relative  importance  of  the  terms,  each  term  has  been 
normalized  for  each  time  step  by  division  of  the  magnitude  of  the  largest 
term  for  that  time  step.  These  four  curves  do  not  show  true  time  histo- 
ries. The  convective  acceleration  term  is  never  dominant. 


Figure  10  is  a plot  of  predicted  and  observed  wave  amplification 
factors  for  a range  of  forcing  wave  periods.  Wave  amplification  is  de- 
fined as  the  bay  level  fluctuation  amplitude,  a^,,  divided  by  the  forcing 
amplitude,  a^.  The  observed  amplification  factors  are  usually  slightly 
smaller  than  the  predicted  values  on  the  Great  Lakes  because  the  temporal 
variation  in  the  prototype  wave  amplitude  produces  transient  effects  that 
do  not  allow  amplification  to  fully  develop. 


The  maximum  amplification  at  Pentwater  occurs  at  a forcing  period  of 
1.8  hours  (Fig.  10)  and  the  maximum  velocity  occurs  at  a period  of  1.4 
hours,  assuming  a constant  forcing  asplitude  of  3 centimeters.  The 
relatively  high  amplification  of  the  forcing  wave,  the  short  Helmholtc 
period,  and  the  large  ratio  between  surface  area  of  the  bay  and  inlet 
cross-sectional  area  (lO**)  cause  relatively  large  inlet  velocities  at 
Pentwater  (0.6  meter  (2  feet)  per  second). 


b.  Tsun^i  Effects  in  a Planned  Inlet.  This  model  may  also  be  used 
to  predict  the  effects  of  a tsunami  at  some  inlets.  This  example  analyzes, 
the  response  of  a hypothetical  inlet-harbor  system  to  an  assumed  tsunami. 
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Response  to  long  wave  excitation  at  Pentwater,  Michigan 
(wave  amplitude  ^ 0.1  foot). 


(1)  Geometry.  The  inlet  is  122  meters  long,  and  24  meters 
(80  feet)  wide  with  a design  depth  of  7.3  meters  (24  feet).  The  bay 
area  is  4.6  x 10®  square  meters  (5  * 10®  square  feet). 

(2)  Forcing.  The  assumed  tsunami  seawater  level  fluctuation  is 
shown  in  Figure  11.  This  tsunami  event  includes  four  waves  with  a period 
of  approximately  20  minutes.  The  third  wave  is  the  largest  with  a crest 
elevation  of  0.4  meter  (1.3  feet). 

(3)  Calibration.  No  data  are  available  for  calibration.  How- 
ever, other  examples  have  shown  that  the  Masch,  Brandes,  and  Reagan 
(1977)  relation  adequately  describes  Manning's  n for  inlets  (eq.  18). 
These  examples  also  show  that  the  empirical  entrance-  and  exit-loss 
coefficients  for  flood  and  ebb  are  approximately  1.0. 

(4)  Results.  The  numerical  model  predicts  that  the  bay  water 
level  fluctuations  will  lag  the  sea  fluctuations  by  approximately  5 
minutes,  and  will  be  smaller  than  the  sea  fluctuations  (Fig.  12).  The 
maximum  flood  and  ebb  velocities  will  be  approximately  2.7  meters  (9  feet) 
per  second  and  will  occur  during  the  third  wave  (Fig.  12).  Because  of 
temporal  acceleration  effects,  the  harbor  will  continue  to  oscillate 
after  the  sea  level  forcing  has  stopped  acting  on  the  harbor.  These 
oscillations  will  die  out  exponentially  due  to  friction  losses  (Fig.  12). 

c.  Masonboro  Inlet,  North  Carolina.  Masonboro  Inlet  is  an  example 
of  a tidal  inlet  connecting  a poorly  defined  bay  to  the  Atlantic  Ocean. 

The  inlet  is  controlled  by  a single  weir  jetty  on  the  northeast  side  of 
the  inlet  (Fig.  13).  The  model  was  calibrated,  using  sea  and  bay  water 
levels  and  inlet  velocities  measured  12  and  13  September  1969. 

(1)  Geometry.  Hydrography  for  1969  was  modeled,  using  a grid 
system  of  four  channels  and  seven  cross  sections  drawn  to  estimate  a flow 
net  for  the  inlet  (Fig.  14).  Grid  measurements  of  cross-sectional  area, 
width,  and  length  were  made  for  each  grid.  Cross-sectional  areas  were 
referenced  to  mean  sea.  level  (MSL),  which  is  1.88  feet  above  the  local 
Beaufort  mean  low  water  datum. 

Although  the  Masonboro  bay  surface  area  is  poorly  defined,  the  tidal 
prism  was  measured.  Therefore,  the  effective  bay  area  is  taken  as  the 
measured  tidal  prism  divided  by  the  mean  bay  tidal  range  measured  at  three 
gages  (gages  3,  4,  and  5 in  Fig.  14)  to  give  an  effective  bay  surface  area 
of  1.8  * 10®  square  feet  at  MSL.  The  bay  surface  area  at  any  point  in 
time,  ^bayt  varies  as  a function  of  the  bay  water  level  because  exten- 
sive lowlana  areas  are  flooded  at  high  water.  The  assumed  relation 
between  bay  surface  area  and  bay  level  is 

■ Aj,(l  ♦ Ph^)  (19) 


with  B - 0.2  estimated  from  hydrographic  charts. 
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Figure  12.  Predicted  inlet-bay  system  response  to  tsunami-generated 
seawater  level  fluctuations. 


The  slopes  of  Masonboro  Inlet  are  estimated  to  be  0.0133  from  the 
hydrographic  survey. 

(2)  Forcing.  Tidewater  levels  measured  outside  the  inlet  at 
0.5-hour  intervals  for  a period  of  25  hours,  12  and  13  September  1969, 
were  used  as  the  input  to  the  model . 

(3)  Calibration.  The  prototype  inlet  velocities  measured  12 

September  1969  were  used  to  calibrate  the  Masonboro  model.  Values  of 
Manning's  n throughout  the  inlet  grid  system  are  assumed  to  be  a func- 
tion of  the  water  depth  below  datum  in  each  grid,  using  the  rela- 

tion recommended  by  Masch,  Brandes,  and  Reagan  (1977): 


n = 0.0377  - 0.000667  (20) 


Flow  is  distributed  throughout  the  inlet  so  that  friction  is  minimized 
in  each  channel  (Hirr*2) . 

Figure  15  shows  the  forcing  tide,  predicted  bay  levels,  inlet  velocity 
at  the  throat,  and  discharge.  The  plot  of  the  importance  of  terms  in  the 
equation  of  motion  shows  that  the  Masonboro  Inlet  flow  is  predominantly 
controlled  by  head  and  frictional  effects.  Temporal  acceleration  becomes 
important  only  for  a short  time  near  slack  water. 

Comparisons  between  observed  and  calibrated  inlet  velocities  show 
that  this  is  a good  model  for  predicting  mean  velocity  in  the  throat  at 
Masonboro  (Fig.  16).  The  standard  deviation  between  observed  and  predicted 
is  0.62. 

This  model  was  also  used  to  estimate  flow  distribution  throughout  the 
inlet  using  output  from  the  subroutine  TABLE  in  the  computer  program. 
Predicted  maximum  ebb  and  flood  velocities  are  illustrated  in  Figures 
17  and  18.  This  calibration  shows  that  highest  velocities  occur  in  the 
deep  gorge  of  the  inlet  adjacent  to  the  jetty. 

d.  Indian  Indian  River  is  an  example  of  a 

single-jettied  inlet  which  connects  an  open  bay  to  the  Atlantic  Ocean 
(Fig.  19).  This  example  illustrates  the  effects  of  a storm  surge  on 
water  motions  in  the  inlet  and  bay. 

(1)  Geometry.  The  inlet,  approximately  1,650  meters  (5,400 
feet)  long,  is  comparatively  straight  and  uniform  throughout  its  length, 
so  the  inlet  is  modeled  using  one  channel  and  five  cross  sections  (four 
grids  long).  The  channel  dimensions  of  1943  (Keulegan,  1967)  were  used 
to  model  the  inlet  geometries  (Fig.  20).  The  bay  area  is  taken  as 
4.2  X 10®  square  feet  (Keulegan,  1967)  and  as  a first  approximation  the 
effect  of  "The  Ditches"  is  neglected  because  flow  through  "The  Ditches" 
is  small  compared  to  the  inlet  discharge  (Keulegan,  1967). 
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Figure  17.  Ebb  velocity  distribution  (minimim  friction  weighting 
1100  hours,  e.s.t.) . 


Figure  18.  Flood  velocity  distribution  (minimum  friction  weighting, 
1800  hours,  e.s.t.). 
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Figure  20.  Schematic  of  Indian  River  model  geometry 


(2)  Forcing  Hydraulics.  The  storm  surge  measured  at  Atlantic 
City,  New  Jersey,  14  and  15  September  1944  is  assumed  representative  of 
the  surge  at  Indian  River  (Jelesnianski,  1967).  This  surge  is  added  to 
the  astronomical  tide  to  obtain  a time  history  of  water  levels  in  the 
Atlantic  Ocean  at  Indian  River  (Fig.  21). 

(3)  Calibration.  No  data  are  available  for  calibration  of  this 
model  to  storm  surge  conditions.  Keulegan  (1967)  found  that  a Manning's 
n of  0.046  applied  to  Indian  River  Inlet,  so  this  value  is  used  in  the 
model . 


(4)  Results.  As  a first  approximation,  this  model  predicts  that 
the  storm  surge  and  astronomical  tide  together  produce  a maximum  flood 
velocity  of  200  centimeters  (6.5  feet)  per  second  followed  by  a maximum 
ebb  velocity  of  110  centimeters  (3.6  feet)  per  second  (Fig.  21).  Bottom 
stress  and  head  are  the  most  important  terms  in  the  equation  of  motion. 
Temporal  acceleration  becomes  important  only  near  slack  water.  Convec- 
tive acceleration  is  negligible.  In  spite  of  the  large  surge  in  the 
Atlantic  Ocean  and  the  high  inlet  velocities,  the  bay  level  changes  are 
relatively  small  because  the  bay  surface  area  is  large  compared  to  the 
inlet  cross-sectional  area  (a  ratio  of  10^)  and  the  surge  only  lasts  for 
a few  hours. 

e.  Cabin  Point  Creek,  Virginia.  Cabin  Point  Creek  is  a shallow, 
natural  tidal  inlet  (Fig.  22,  inlet  1)  connected  to  the  Potomac  River. 

This  example  shows  how  the  model  can  be  used  to  assist  in  the  design  of 
a second  inlet  (inlet  2). 

(1)  Geometry.  The  area  of  the  bay,  measured  from  hydrographic 
charts,  is  3.3  x 10®  square  meters  (3.5  x 10®  square  feet)  and  the  length 
of  the  natural  inlet  is  580  meters  (1,900  feet).  Only  one  cross  section, 
located  one-third  of  the  length  of  the  inlet  from  the  sea,  was  measured 
for  the  inlet  in  May  1976  (Seelig,  1976).  The  maximum  water  depth  was 
0.6  meter  (2.1  feet)  with  an  inlet  width  of  15  meters  (50  feet).  This 
inlet  was  modeled,  using  a three-channel  flow  net.  The  bay  surface  varia- 
tion parameter,  B,  was  determined  from  hydrographic  charts  as  0.2. 

(2)  Forcing.  Water  level  measurements  in  the  Potomac  River  for 
an  18-hour  interval  on  24  and  25  May  1976  were  used  to  force  the  model. 

(3)  Calibration.  The  computer  program  model  of  Cabin  Point  was 
calibrated  by  comparing  observed  and  predicted  bay  levels.  Masch,  Brandes, 
and  Reagan's  (1977)  recommended  relation  between  n and  water  depth  for 
depths  less  than  1.2  meters  and  ebb-  and  flood- loss  coefficients  of  1.0  to 
adequately  model  bay  levels  (Pig.  23). 

(4)  Prediction.  The  calibrated  model  was  used  to  predict  the 
hydraulics  of  the  system  if  a second  Inlet  (hypothetical),  91  meters 
(300  feet)  long,  1.2  meters  deep,  and  15  meters  wide,  were  cut  at 
location  2 shown  in  Figure  22.  The  model  predicts  that  the  bay  tide 
range  would  increase  by  a factor  of  four  (from  0.36  to  1.49  feet,  0.1  to 
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Figure  22.  Cabin  Point  Creek,  Virginia 


0.45  meter)  and  the  velocities  in  inlet  1 would  decrease  (Table  1). 
The  maximum  water  velocities  in  inlet  2 would  be  -1.3  and  1.7  feet 
(-0.4  and  0.5  meter)  per  second  for  ebb  and  flood. 


Table.  Predicted  Cabin  Point  Creek  hydraulics.^ 


Tide 

range 

Max . ebb 

Max.  flood 

Cabin  Point  Creek 

bay  1 river 
(ft) 

velocity 

(ft/s) 

velocity 

(ft/s) 

Inlet  1 

(24  and  25  May  1976) 

0.36 

1.49 

-0.6 

0.9 

Predicted  with  second  inl^t 
Inlet  1 1 

Inlet  2^  1 

\ 

-1.49 

-1.49 

1.49 

1.49 

-0.3 

-1.3 

0.3 

1.7 

^ Inlet  2 has  L = 300  feet,  B = 50  feet,  and  D = 4 feet. 


To  illustrate  the  effect  of  the  width  of  inlet  2 on  velocities  and 
bay  levels,  the  model  was  run  with  several  widths  between  0.6  meter  and 
30  meters  (100  feet);  the  results  are  plotted  in  Figure  24.  This  analysis 
shows  that  as  the  width  of  inlet  2 is  increased  the  bay  tide  range  in- 
creases as  the  bay  tide  more  closely  follows  the  river  tide.  At  a width 
of  12  meters  (40  feet),  the  bay  range  is  the  same  as  the  river  range. 

As  the  width  of  inlet  2 increases  and  the  bay  tide  more  closely  approxi- 
mates the  sea  tide,  tjie  head  differentials  along  the  inlets  decrease, 
so  the  inlet  velocities  decrease  (Fig.  24).  The  plot  shows  that  for  a 
9-meter-wide,  small  Boat  navigation  channel  at  inlet  2,  the  inlet  veloc- 
ities will  be  small  *(^ess  than  or  equal  to  2 feet  per  second)  for  the 
given  tide  condition.  |,See  Seelig  (1976)  for  additional  details  of  the 
Cabin  Point  Creek  example. 

V.  SUMMARY  AND  CONCLUSIONS 

A numerical  model  based  on  an  area  averaged  momentum  equation  for 
the  inlet,  and  a continuity  equation  for  the  bay  has  been  shown  to  give 
good  predictions  of  bay  levels  and  mean  inlet  velocities  for  a variety 
of  inlets  and  forcing  conditions.  The  model  is  designed  for  cases  where 
the  bay  water  level  fluctuates  uniformly  throughout  the  bay  and  the 
volume  of  water  stored  in  the  inlet  between  high  and  low  water  is  negli- 
gible compared  to  the  prism  of  water  that  moves  through  the  inlet. 

The  ease  of  use  and  low  cost  of  this  model  make  it  ideal  for  obtaining 
a first  estimate  of  inlet  hydraulics  for  inlets  forced  by  the  astronomical 
tide,  storm  surge,  lake  seiching,  and  tsunamis.  The  computer  program 
INLET  (App.  D) , based  on  the  model,  can  be  used  for  hydraulic  calculations 
when  one  or  more  inlets  connect  a bay  to  a sea. 
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Figure  24.  Effect  of  the  width  of  inlet  2 on  inlet  and  bay  hydraulics 
Cabin  Point  Creek,  Virginia. 
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There  are  several  ways  that  this  program  can  be  used  to  aid  the  design 
engineer  and  researcher.  For  example,  when  designing  an  inlet  system 
many  physical  parameters  can  be  easily  varied  in  the  computer  program  to 
determine  the  influence  of  each  factor  on  inlet  and  bay  hydraulics. 

The  model  can  also  be  used  to  obtain  an  estimate  of  the  magnitude  of 
measurements  to  be  made  in  the  field.  For  example,  at  Masonboro  the  model 
indicates  that  small  head  differences  generate  large  inlet  velocities. 

This  means  that  Masonboro  water  level  measurements  should  be  exceptionally 
accurate  and  precise.  The  model  shows  that  for  Pentwater,  waves  with 
periods  of  1 to  2 hours  and  amplitudes  of  0.1  foot  can  generate  signifi- 
cant inlet  velocities.  This  suggests  that  the  water  level  recorder  and 
stilling  well  (Seelig,  1977)  at  Pentwater  should  be  designed  to  measure 
the  long  waves  with  precision. 

This  numerical  model  can  be  used  to  examine  the  relative  importance 
of  terms  in  the  equation  of  motion  of  inlet  hydraulics.  For  example, 
application  of  the  model  for  several  locations  shows  that  at  tidal  inlets, 
pressure  head  and  friction  are  the  dominant  terms  in  the  equation  of 
motion  throughout  most  of  the  tide  cycle.  Temporal  acceleration  becomes 
important  only  near  slack  water.  At  nontidal  inlets,  such  as  on  the 
Great  Lakes,  three  terms  (temporal  acceleration,  pressure  head,  and 
friction)  are  important  throughout  a major  part  of  the  water  level  cycle. 
The  strong  influence  of  temporal  acceleration  on  Great  Lakes  inlets 
causes  amplification  of  long  waves  near  the  natural  pumping  period  of  the 
inlet-bay  system. 

In  addition  to  INLET,  a second  computer  program,  INLET2,  is  now 
available.  INLET2  can  be  used  to  predict  hydraulics  for  systems  consist- 
ing of  one  or  more  inlets  bays,  and  seas,  including  complex  systems. 

These  computer  decks  may  be  obtained  from  the  Automatic  Data  Processing 
Coordinator  (CERDP),  Coastal  Engineering  Research  Center,  Fort  Belvoir, 
Virginia  22060. 

VI.  RECOMMENDATIONS 

There  are  several  changes  recommended  for  the  computer  programs 
INLET  and  INLET2  to  increase  the  generality  of  the  model.  Storage  of 
water  in  inlets,  and  radiation  and  wind  stresses  may  be  added  to  the 
equations  of  motion  for  the  inlet.  Other  weighting  functions  for  dis- 
tributing flow,  for  example  a weighting  function  which  includes  the  jet 
action  of  water  exiting  the  inlet,  may  be  added  to  the  program.  Higher 
order  relations  to  evaluate  bay  surface  area,  inlet  cross-sectional  area, 
and  inflow  into  the  bay  from  sources  other  than  the  inlet  can  also  be 
added  to  the  computer  program.  A bottom-friction  relation  designed  for 
unsteady  flow  may  slightly  increase  the  quality  of  the  results. 
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APPENDIX  A 


WEIGHTING  FUNCTIONS  TO  MINIMIZE  FRICTION 


Let  the  total  bottom  stress  in  an  inlet  cross  section  be  described  by 


where  C^  is  a parameter  that  can  be  determined  using  results  from  the 
previous  time  step: 


where  W^  is  an  unknown  weighting  function  and  Q is  the  inlet  discharge 
C^  can  be  evaluated  for  each  grid  (as  in  the  computer  subroutine  WTl),  or 
for  each  channel  (WT2).  In  many  cases,  it  is  reasonable  to  assume  that 
flow  in  an  inlet  is  distributed  so  that  total  friction  is  minimized  (i.e., 
fiOW  will  follow  the  path  of  least  resistance).  In  that  case,  find 
to  minimize  F^.  To  find  W^  set  up  a new  function,  F^,  which  includes 
a Lagrangian  multiplier.  A: 


Then  differentiate  with  respect  to  each  W^  to  obtain  a series  of 
equations: 


Also .differentiate  equation  (A-3)  with  respect  to  X to  obtain 


Setting  equations  (A-4)  and  (A-5)  to  lero  will  give  the  values  of 
to  obtain  the  minimum  or 
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(A- 6) 


Substituting  equation  (A-6)  into  equation  (A-5)  for  W^: 


1 = 0 


and  solving  for  A: 


A = 


(A-7) 


CA-8) 


Then  substitute  X into  equation  (A-6)  to  obtain  the  weighting  function 
to  minimize  friction: 


(A-9) 


where  is  given  by  equation  (A-2). 

This  function  to  determine  the  weighting  for  minimum  friction  can 
be  applied  in  several  ways.  If  flow  is  allowed  to  cross  channels,  the 
minimum  friction  weighting  can  be  applied  at  each  cross  section  of  grids 
(WTl).  This  will  give  the  lowest  possible  friction  in  the  aodel.  If 
flow  is  not  allowed  to  cross  channels,  then  the  routine  can  be  applied 
to  distribute  flow  in  each  channel  to  minimize  friction  (NT2).  This 
will  give  more  frictional  resistance  than  in  NTl. 


APPENDIX  B 


CONSTRUCTING  FLOW  NETS 


A flow  net  for  an  inlet  is  a series  of  subchannels  and  cross  sections 
that  divide  the  inlet  into  a set  of  grid  cells.  The  primary  purpose  of 
the  flow  net  is  to  evaluate  bottom  friction  throughout  the  inlet  at  each 
time  step.  The  process  of  constructing  the  flow  net  is  subjective.  The 
subjectivity  can  be  reduced  by  following  the  procedures  described  here. 
Friction  in  each  cell  is  determined  by  taking  the  mean  geometric  and 
hydraulic  conditions  throughout  each  cell  and  applying  them  at  the  cen- 
troid of  the  cell  through  the  use  of  Manning's  bottom-friction  relation- 
ship. Total  inlet  friction  is  the  sum  of  friction  in  all  the  cells. 

The  flow  net  should  be  drawn  so  that  flow  parallels  the  subchannels  and 
is  perpendicular  to  the  cross  sections.  In  prototype  inlets,  the  actual 
flow  net  system  changes  with  time  and  these  changes  are  partially  accounted 
for  in  the  model  by  reevaluating  the  weighting  function  (IWT*1.  and  Ikrr=2) 
at  each  time  step. 

1.  Drawing  Flow  Nets. 

One  method  of  drawing  a flow  net  for  an  inlet  is  by  assuming  that 
flow  is  distributed  at  each  cross  section  so  that  friction  is  minimized 
(1WT=1  in  the  computer  program;  see  App.  A).  The  first  procedure  for 
applying  this  method  is  to  (a)  draw  a number  of  cross  sections  which  are 
approximately  perpendicular  to  flow  for  the  inlet,  (b)  apply  the  minimum 
friction  weighting  function  to  each  section  to  determine  channel  locations 
so  that  flow  through  each  channel  will  be  equal,  and  (c)  check  to  see  that 
the  estimated  cross  sections  are  perpendicular  to  the  channels  drawn  in 
(b).  If  not,  repeat  steps  (a)  and  (b) . 

A hypothetical  inlet  is  used  to  illustrate  the  technique  for  drawing 
flow  nets  (Fig.  B-1).  This  idealized  inlet  has  a maximum  depth  of  5.5 
meters  (18  feet)  below  MSL  in  the  throat,  is  366  meters  (1,200  feet)  wide 
at  the  minimum  width,  and  is  symmetric  about  the  longitudinal  axis  of  the 
inlet.  A deep  gorge  runs  through  the  centerline  of  the  inlet. 

The  first  step  in  drawing  a flow  net  grid  system  is  to  draw  a number 
of  inlet  cross  sections  which  are  approximately  perpendicular  to  flow  in 
the  inlet.  These  estimated  cross  sections  are  drawn,  keeping  in  mind 
that  in  the  throat  and  near  the  centerline  of  the  inlet  most  flow  will  be 
parallel  to  the  main  axis  of  the  inlet.  Beyond  the  inlet,  the  flow  will 
move  in  more  of  a radial  direction  where  the  inlet  has  the  effect  of  a 
source.  Figure  B-2  illustrates  15  cross  sections  estimated  for  the  ideal- 
ized inlet.  These  sections  are  placed  closer  together  in  the  throat  of 
the  inlet  because  this  is  a high  friction  region. 

The  second  step  is  to  determine  water  depths  at  points  along  each 
cross  section,  and  based  on  this  information,  determine  channel  loca- 
tions so  that  flow  is  distributed  to  minimize  the  head  loss  due  to 
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friction.  This  task  is  aided  by  the  FORTRAN  computer  program,  NET, 
discussed  later  in  this  appendix.  Figure  B-3  shows  the  channel  locations 
estimated  so  that  flow  is  equal  in  all  channels  for  section  12  of  the 
idealized  inlet.  Channels  are  widely  spaced  in  shallow  areas  of  the 
inlet  because  friction  in  these  areas  is  high  and  the  flow  per  unit  cross- 
section  width  is  small.  Channels  have  much  narrower  spacing  in  deep  areas 
because  friction  is  low  and  the  greater  depth  allows  for  more  flow  per 
unit  width. 

Channel  location  computations  are  made  for  all  cross  sections  and  the 
locus  of  channel  locations  drawn  to  form  a flow  net  (Fig.  B-4). 

Another  output  from  the  program  NET  is  the  predicted  relative  amount 
of  frictional  dissipation  per  unit  inlet  length  for  each  of  the  cross 
sections.  This  friction  loss  is  plotted  on  the  left  side  of  Figure  B-4. 
The  flow  net  evaluates  friction,  so  those  sections  located  far  from  the 
inlet  where  friction  is  small  can  be  eliminated.  For  example,  at  sections 
1 and  15,  friction  is  extremely  small  because  the  water  velocities  at 
these  points  are  low.  As  a rule-of-thumb  to  account  for  most  friction 
losses,  the  first  and  last  cross  sections  of  a flow  net  should  have  cross- 
sectional  areas  approximately  10  times  larger  than  the  minimum  area  in 
the  throat. 

2.  Documentation  of  the  Computer  Program  NET. 

NET  is  a FORTRAN  computer  program  to  aid  the  user  in  drawing  a flow 
net  grid  system  of  an  inlet.  Table  B-1  lists  the  information  needed  for 
input  to  the  program.  This  information,  including  cross-sectional  geome- 
try, is  analyzed  using  the  weighting  function  to  minimize  friction  (App. 

A)  to  determine  the  location  of  channel  boundaries  so  that  flow  through 
all  channels  will  be  equal. 

Output  of  the  program  includes  a printer  plot  of  each  cross  section 
and  the  relative  amount  of  discharge  at  points  along  the  section  (Fig. 
B-3).  Dimensions  of  channels  for  equal  flow  in  all  channels  are  plotted 
and  output  in  tabular  form  to  be  used  as  input  to  INLET  (Fig.  B-3).  A 
summary  table  gives  the  relative  friction  loss  per  foot  of  channel  length 
for  each  cross  section  (Table  B-2).  This  table  shows  which  areas  of  the 
inlet  are  most  important  in  dissipating  friction;  e.g.,  section  8 has  the 
highest  friction  loss  for  the  idealized  inlet  (Table  B-2,  Fig.  B-4). 

Table  B-3  is  a listing  of  the  program  NET  including  comments  on 
program  flow. 

Output  information  from  NET  showing  channel  location  (e.g..  Fig.  B-3) 
should  be  added  to  each  cross  section  and  channels  drawn  by  connecting 
the  channel  location  points  (e.g.,  points  A,  B,  C,  D) . If  the  cross  sec- 
tions are  not  perpendicular  to  the  computed  channel  locations  in  high 
friction  areas  of  the  inlet,  these  cross  sections  can  be  redrawn  and 
channel  locations  recomputed  to  obtain  a more  accurate  flow  net.  How- 
ever, this  should  not  be  necessary  in  most  cases.  Point  A in  Figure  B-4 
shows  that  section  1 is  not  perpendicular  to  flow  and  it  should  be  re- 
drawn if  used  in  a model. 
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Figure  B-3.  Computed  channel  locations  for  a cross  section  using 
minimum  friction  weighting  (output  from  program  NET) . 
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Figure  B-4.  Idealized  inlet  flow  net  and  predicted  friction  loss 
along  the  inlet  length  (adapted  from  output  from  the 
computer  program  NET). 


Table  B-1.  Input  to  the  computer  program  NET. 


Card 

Format 

Description 

1 

no 

Read  the  number  of  inlets  to  be  input. 

2 

no,  FIO.O 

Read  number  of  channels,  inlet  discharge 
(cubic  feet  per  second) . 

3 

no,  FIO.O,  A4 

Read  number  of  depth  readings,  spacing  between 
depth  readings  (feet),  title  of  section. 

Spacing  of  depth  readings  should  be  small 
enough  so  that  linear  interpolation  between 
readings  will  adequately  describe  the  bottom 
topography  for  each  segment. 

4 

16F5.1 

Read  depths  at  equal  spacing  across  the  channel 
(feet).  The  first  and  last  values  should  be 

0.0  (at  the  waterline  at  each  end  of  the  sec- 
tion). Card  type  4 is  repeated  as  many  times 
as  necessary  to  read  all  the  depth  values  for 
this  section  (e.g.,  for  20  depth  readings 
there  would  be  two  card  type  4's  for  that 
section) . 

NOTE. --For  each  cross  section  for  an  inlet  repeat  card  types  3 and  4. 


Place  a blank  card  after  the  last  cross  section  of  an  inlet  to  in- 
dicate the  end  of  that  inlet. 

For  each  inlet  repeat  card  types  2,  3,  and  4. 


Table  B-2.  Suamary  of  friction  losses. 


Section 

Loss/ ft  channel  length 
(dimensionless) 

1 

0.206 

2 

0.398 

3 

0.699 

4 

1.84S 

5 

5.740 

6 

11.298 

7 

13.057 

8 

21.451 

9 

14.598 

10 

14.222 

11 

6.860 

12 

4.360 

13 

2.471 

14 

1.695 

IS 

1.100 
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Table  B-3.  Listing  of  the  computer  program  NET. 
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A suggested  flow  net  grid  system  which  includes  important  friction 
zones  of  the  idealized  inlet  is  shown  in  Figure  B-5. 
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Figure  B-S.  A suggested  grid  cell  system  for  the  idealized  inlet 
(seven  channels  and  eight  cross  sections). 


APPENDIX  C 


ALTERNATIVE  DEVELOPMENT  OF  THE  CONVECTIVE  ACCELERATION  TERM 

In  equation  (1)  it  was  assumed  that  the  convective  acceleration  term 
could  be  integrated  along  the  length  of  the  inlet.  The  resulting  expres- 
sion for  convective  acceleration  is  1/2  (l/A?  - 1/A|)  (see  eq.  6,  this 
1*  option  ICONV  *0  in  the  computer  program  INLET).  This  derivation 
assumes  that  the  effective  head  converted  into  acceleration  of  water  as 
water  enters  the  inlet  can  be  regained  as  the  water  exits  the  inlet. 
However,  some  of  this  head  may  be  lost  due  to  turbulence  and  eddies. 

If  it  is  assumed  that  convective  acceleration  is  described  by  a 
coad>ined  empirical  exit-  and  entrance-loss  coefficient,  Cj)  (Keulegan, 
1967),  and  the  water  velocity  far  from  the  inlet  is  small,  tnen  the  equa- 
tion of  motion  for  the  inlet  (eq.  12)  can  be  rewritten: 


dQ  "'s 

* "7  /Wn' 


CC-1) 


where  A,n^„  is  the  minimum  inlet  cross-sectional  area  (the  area  at  the 
throat).  If  an  inlet  model  corresponding  to  equation  (C-1)  is  desired, 
the  option  ICONV  ^ 0 should  be  used  in  the  computer  program  INLET. 

The  values  of  Cp  for  both  ebb  and  flood  must  also  be  specified 
(CDE  and  CDF  input  to  INLET).  As  a first  approximation,  values  of  Cq 
should  be  about  1.0  for  ebb  and  flood  for  rapidly  converging  and  diverg- 
ing inlets.  This  means  that  one  velocity  head  is  lost  due  to  the  con- 
traction and  expansion  process  through  the  inlet.  However,  if  the  inlet 
has  a different  shape  on  the  sea  and  bay  ends,  this  loss  coefficient  may 
be  different  for  ebb  and  flood.  The  hypothetical  inlet  in  Figure  C-1 
illustrates  possible  reasons  for  different  loss  coefficients  for  ebb  and 
flood.  On  floodtide,  the  water  converges  in  a radial  fashion  and  acceler- 
ates as  it  approaches  the  throat  of  the  inlet.  Convective  acceleration 
develops  one  throat  velocity  head  at  the  throat.  As  the  water  exits  the 
inlet  a significant  part  of  this  velocity  head  is  lost  in  the  form  of 
eddies.  The  loss  coefficient,  Cjj,  for  floodflow  is  expected  to  be 
equal  or  exceed  1.0  for  this  exaa^le.  ^ 

On  the  ebbtide,  water  approaches  the  inlet  in  an  approximately  radial 
fashion  and  then  diverges  slowly  due  to  the  funnel  shape  of  the  seaside 
of  the  inlet.  This  pattern  of  flow  causes  little  head  loss  due  to  eddies, 
so  a significant  part  of  head  converted  into  acceleration  of  water  on 
entrance  is  regained  as  the  flow  exits  the  inlet.  The  value  of  Cp  for 
this  flow  condition  is  e]q)ected  to  be  less  than  1.0. 

As  a first  approximation  values  of  Cp  for  ebb  and  flood  can  be 
taken  as  1.0  for  rapidly  converging  inlets.  These  loss  coefficients  can 
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Figure  C-1.  Flood  and  ebb  convective  losses  for  a hypothetical  inlet 
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then  be  adjusted  during  calibration  to  account  for  differences  in  ebb 
flow  and  floodflow.  For  example,  if  in  the  prototype  flood  velocities 
are  lower  than  ebb  velocities,  the  values  of  the  ebb-  and  flood-loss 
coefficients  should  be  adjusted  until  the  model  reproduces  observed 
velocities.  A larger  value  of  the  loss  coefficient  will  mean  a greater 
head  loss  and  associated  lower  velocities. 


APPENDIX  D 


COMPUTER  PROGRAM  (INLET)  DOCUMENTATION 


1 . Program  Description. 

The  numerical  model  to  predict  inlet  hydraulics  described  in  the 
text  is  programed  in  FORTRAN  for  a CDC  6600  computer.  The  simultaneous 
differential  equations  are  solved  using  a variable  time  step  Runge-Kutta- 
Gill  marching  procedure.  The  organization  of  the  computer  program  is 
shown  in  Figure  D-1.  The  following  is  a brief  description  of  each  routine: 

INLET  is  the  main  routine  which  controls  input -output  and  calls  sub- 
routines  to  execute  a specific  task.  Figure  D-1  sunAiarizes  control 
throughout  the  program.  Variable  names  in  the  program  corresponding  to 
symbols  used  in  the  text  are  given  in  Symbols  and  Definitions.  The 
program  is  organized  to  accept  up  to  three  inlets  connecting  the  bay  to 
the  sea,  up  to  seven  channels  for  each  inlet,  and  up  to  eight  cross 
sections  (seven  grids  long). 

Subroutine  HELM  uses  an  iterative  method  of  estimating  the  natural 
pumping  period  or  Helmholtz  period,  Tg',  for  the  inlet-bay  system  by 
neglecting  friction  in  the  inlet  to  give 


T • - 2,  J35ZE1S  , 


where  L'  is  added  inlet  length  due  to  radiation,  and  given  by: 

L*  - ^ In  LllA  (D-2) 

• (TJSTij 

Subroutine  RKCS  is  a routine  to  solve  simultaneous  differential 
equations.  This  subroutine  was  adapted  from  the  Scientific  Subroutine 
Package  (International  Business  Machines,  1970). 

Subroutine  SETEQ  evaluates  the  right-hand  side  of  the  equation  of 
motion,  one  for  each  inlet,  and  the  continuity  equation  between  the  inlet 
and  bay  for  each  time  step.  This  routine  also  evaluates  the  relative 
rank  of  the  four  terms  in  the  equation  of  motion  for  flow  in  each  inlet. 

Subroutine  LEVEL  determines  the  water  level  in  the  grids  at  each  time 
step.  The  routine  interpolates  the  level  between  the  sea  and  bay  based 
on  the  relative  amount  of  friction  in  each  grid  cell. 
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Subroutine  TPWRTE  writes  hydraulic  results  from  each  time  step  on  a 
tape  or  disc,  so  that  this  information  can  be  used  later  by  the  output 
routines. 

Subroutine  TABLE  outputs  a table  of  instantaneous  hydraulics  each 
time  the  routine  is  called. 

Subroutine  SEA  determines  the  water  level  in  the  sea  as  a function  of 
time  either  for  a given  sine  wave  or  by  interpolating  equal  time-series 
data. 

Subroutine  WTl  determines  the  grid  weighting  function  by  assuming 
that  the  flow  is  distributed  across  each  section  so  that  friction  is 
minimized.  This  routine  allows  flow  to  cross  channel  boundaries,  but 
assumes  that  this  flow  will  be  small,  so  the  flow  is  neglected  in  the 
equation  of  motion.  See  Appendix  A for  a derivation  of  the  procedure 
used  in  this  routine. 


Subroutine  WT2  is  similar  to  WTl,  except  that  flow  is  not  allowed  to 
cross  channel  boundaries,  and  the  flow  is  distributed  in  each  channel  so 
that  friction  is  minimized. 

Subroutine  WT3  determines  the  weighting  function  so  that  flow  is 
distributed  equally  in  all  grids.  This  is  generally  unrealistic,  since 
it  will  be  difficult  to  visually  draw  this  grid  system.  However,  this 
routine  is  useful  since  it  provides  an  upper  limit  on  frictional  effects 
and  therefore  gives  a lower  limit  of  bay  levels  and  inlet  velocities. 

This  weighting  can  be  used  to  model  simple  geometry  inlets  where  only  one 
channel  is  used  to  represent  the  inlet. 

Subroutine  CRIT  prints  a table  of  critical  instantaneous  hydraulics 
(i.e.,  at  high  water,  low  water,  maximum  velocity,  and  maximum  discharge). 
This  table  is  determined  by  storing  a summary  of  conditions  for  each  time 
step,  then  scanning  this  list  for  critical  values  (routine  adapted  from 
Huval  and  Wintergerst,  1977), 

Subroutine  GRPHC  plots  mean  inlet  hydraulics  by  scaling  hydraulics  in 
storage  and  plotting  the  time  interval  requested  on  a digital  x-y  pen 
plotter  (adapted  from  Huval  and  Wintergerst,  1977). 

Subroutine  READ IN  is  used  by  GRPHC  to  read  data  in  storage  and  scale 
values  for  plotting  (adapted  from  Huval  and  Wintergerst,  1977). 


The  computer  program  INLET  requires  the  following  input,  one  deck 
for  each  inlet-bay  system: 
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• V.4' 


r 


Card 

Type  Variables  Format 

1 AUBLl  4A10 

ALABL2  4A10 

2 5110,  2F10.5,  no 
NINLET 

I PLOT 
IWT 

ITABLE 

Cl,  C2 

ICONV 

3 3F10.5,  E10.4, 

3F10.5,  2F5.1 

T 

DELT 

AO 

AB 

BETA 

ZETA 


Description 

First  line  of  title 
Second  line  of  the  title 


Number  of  inlets 
=1  for  plot  of  results 
Weighting  type 

=1  flow  distributed  to  minimize 
friction  at  each  section 

=2  flow  distributed  in  each  chan- 
nel to  minimize  friction 

s3  equal  flow  in  channels 

=1  for  tables  of  instantaneous 
hydraulics 

Manning's  n evaluated  by: 
n=Cl-C2*D;  where  D is 
Stillwater  depth.  If  blank, 
default  values  of  Cl  = 0.03777 
and  C2  = 0.000667  are  assumed. 

ICONV  = 0 for  convective  acceler- 
ation as  derived  in  the  text 
ICONV  f 0 for  alternative 
derivation  (see  App.  C). 


Forcing  period,  hours 

Approximate  time  increment 

Forcing  wave  amplitude  (feet) 

Bay  area  at  datum  (square  feet) 

Bay  area  variation  parameter 

Inlet  side  slope 
D (z)/D(y) 


% 


( 


•f 

■'5 
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Variables 


Format 


Description 


Card 

im. 

QINFLO 

CDF 

CDE 

4 2110,  FIO.O 
IC 

IS 

QINT 

5 (one  card  lOX,  7F10.5 

per  sec- 
tion) 

A' 

6 (one  card  lOX,  7F10.5 

per  sec- 
tion) 

B' 

7 (one  less  lOX,  7F10.5 

card  than 

sections) 

L* 


Bay  inflow  from  sources  other 
than  the  inlet  (cubic  feet  per 
second) 

An  empirical  flood-loss  coeffi- 
cient 

An  empirical  ebb-loss  coefficient 


Number  of  channels 

Number  of  cross  sections 

Estimated  inlet  discharge  at 
the  time  the  model  starts 


Cell  cross-sectional  areas  at 
the  ends  of  each  cell  at  datum 
(square  feet)  (see  Fig.  B-2) 


Grid  cell  widths  for  the  end  of 
each  cell  (feet)  (see  Fig.  2) 


Lengths  of  the  sides  of  cells 
(see  Fig.  2)  (one  less  card 
than  number  of  sections,  one 
more  value  per  card  than  the 
number  of  channels) 


For  card  types  5,  6,  and  7 there  will  be  one  card  for  each  cross 
section  of  the  inlet.  The  first  card  will  be  for  the  first  cross  section, 
i.e.,  the  section  closest  to  the  sea;  the  last  section  is  adjacent  to  the 
bay.  The  first  value  on  each  card  will  correspond  to  the  first  channel 
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Card 

Type  Variables  Format 

which  is  adjacent  to  land,  and  the  last 
last  channel,  which  is  also  adjacent  to 

(FOR  MORE  THAN  ONE  INLET  CONNECTING 
TYPES  3 TO  7 FOR  EACH  ADDITIONAL  INLET) 


8 

TDEL 

34X,  F6.2 

NPTS 

6X,  13 

9 

(optional -no 

cards  if  NPTS 

Y 


10  (optional-two  plot  cards,  first 

on  card  type  1) 

8F10.5/,3F10.5,  110 


Description 

value  on  a card  will  be  for  the 
land  (Fig.  2). 

THE  BAY  TO  THE  SEA  REPEAT  CARD 


Water  level  sampling  interval 
(minutes) 

Number  of  sample  points 
= 0 for  no  data 

from  card  type  8) 

Eight  water  level  values  per  card, 
as  many  cards  to  include  NPTS 
points;  start  the  model  at  a time 
when  the  sea  level  is  zero.  Use 
25  or  more  points  per  forcing 
cycle  for  best  results  (i.e., 
levels  at  30-  or  15-minute  inter- 
vals for  a 12-hour  tide). 

card  used  only  if  IPLOT  = 1 


XO 

Starting  time  of  plot  (hours) 

XF 

Ending  time  of  plot  (hours) 

SCALX 

Time  scale  (hours  per  inch) 

YLO 

Minimum  value  of  water  levels 
(feet) 

YL 

Overall  height  of  plot  (inches) 

YLSCAL 

Scale  of  water  level  height 
(feet  per  inch) 

YRO 

Minimum  flows  (thousands  cubic 
feet  per  second) 

YRSCAL 

Scale  of  flows  (thousands  cubic 
feet  per  second  per  inch) 
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Card 

Jm. 


Variables 


Format 


Description 


Second 

card 


Y\'0 

Minimum  velocity  (feet  per  second) 

YVSCAL 

Scale  of  velocities  (feet  per 

second  per  inch) 

SCALE 

Scale  factor  for  total  plot  size 

IQ 

IQ  = 0 for  no  plot  of  inlet 

discharge 

11  If  a plot  is  requested,  repeat  card  types  8 and  9 for  observed 
bay  levels  to  compare  with  predictions  (card  type  8 required; 
use  NPTS  = 0 for  no  observed  bay  levels) . Only  one  set  of  card 
types  10  and  11  will  be  required  for  plotting  even  though  the 
system  modeled  may  have  more  than  one  inlet. 

12  End  of  file  card. 

The  inlet  data  for  a computer  run  of  Masonboro  Inlet  are  shown  in 
Figure  D-2. 

3.  Output. 

The  types  of  output  include  (a)  a summary  table  of  grid  dimensions, 
input  parameters,  and  the  Helmholtz  period  of  the  system  estimated,  assum- 
ing there  is  no  friction  in  the  inlet;  (b)  (optional)  summary  tables  of 
instantaneous  inlet  hydraulics;  (c)  (optional)  a pen  plot  of  inlet 
hydraulics;  and  (d^  a table  summarizing  critical  points  throughout  model 
operation,  such  as  high  water,  low  water,  point  of  maximum  discharge, 
and  maxim'jm  velocity.  Samples  of  input  and  output  for  the  Masonboro  run 
are  given  in  Figures  D-3,  D-4,  and  D-5.  The  computer  plot  ia  included 
in  the  text  (Fig.  15). 

4 . Computer  Program. 

A listing  of  the  computer  program  INLET  follows  the  sample  output 
flekle  0-S).  The  program  was  written  in  FORTRAN  IV  for  a CDC  6600  com- 
peer Bith  plotter.  Control  cards,  plotting  instructions, and  file  con- 
tmmtp  msf  liave  to  be  changed  for  other  computers.  If  no  plotter  is  avail - 
thtm  am^rowtifie  caiPNC  and  the  call  to  the  subroutine  in  the  main 
gyiypH,  mmr  remmved. 
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Figure  D-2.  Sample  of  input  data  for  a computer  run  of 


Masonboro  Inlet,  North  Carolina. 
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Figure  D-3.  Sanple  output  from  INLET  (summary  table  for 
Masonboro  Inlet  input  data) . 
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Figure  D-4.  Sample  output  from  INLET  (siimmary  table  of  instantaneous 
hydraulics  for  Masonboro  after  6 hours  of  model  time). 
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Figure  D-5.  Sample  output  from  INLET  (table  of  critical  points  for 
the  model  time:  high  water,  low  water,  etc.,  for 
Masonboro  Inlet). 
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Table  D-1,  Listing  of  the  computer  program  INLET 
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